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Abstract 

Starting from a Caldeira-Leggett model for the interaction of a system 
with an environment, Joichi, Matsumoto, and Yoshimura have reconsid- 
ered the derivation of the quantum Boltzmann equation. They find an 
extra term that accounts for the effects of virtual particles, and which 
drastically changes the results for relic densities of stable, weakly inter- 
acting massive particles (WIMPs), and for the decay products of unsta- 
ble particles. We show, however, that this modified Boltzmann equation 
does not properly account for the interaction energy between the mas- 
sive particles (which are decaying or annihilating) and the thermal bath 
of light particles. We argue that the conventional Boltzmann equation 
gives the correct result. 
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Joichi, Matsumoto, and Yoshimura (hereafter JMY) [|l],0] and Matsumoto and Yoshimura 
(hereafter MY) 0] have carefully reconsidered the derivation of the quantum Boltzmann 
equation for heavy particles embedded in a thermal bath of light particles. JMY treat 
the case of unstable heavy particles, and MY treat the case of a stable, weakly interacting 
massive particles (WIMPs) that annihilate to light particles. In both cases they find new 
terms in the quantum Boltzmann equation that account for the effects of virtual heavy par- 
ticles, and that drastically change the usual formula for the equilibrium abundance of these 
particles. In particular, the usual calculation of relic abundances of WIMPs is completely 
changed, with the result that a WIMP with a mass in excess of approximately 1 GeV is would 
overdose the universe, for a broad range of interaction strengths with the light particles. 

These surprising results must be taken seriously, since previous derivations of the quan- 
tum Boltzmann equation (for weakly interacting massive particles) can involve uncontrolled 
approximations and possibly arguable assumptions (see, e.g., 0] for a typical treatment). In 
this context the analyses of JMY and MY are among the most rigorous ones available. 

The basic issue raised by JMY is more easily understood in the context of an unstable, 
decaying particle (rather than stable, annihilating particles). Given a spin- zero particle ip 
with mass M and energy -E(p) = (p^ + M^)^/^ at a temperature T, the conventional formula 
for its equilibrium number density is 



J (27r)3 e^(p)/^ - 1 

_/C(3)TV7r2 for T»M, , . 

\ (MT/27r)3/2e-*VT for T < M. ^ ' 

JMY, on the other hand, argue that if the particle is weakly coupled, and unstable with a 
decay width F <^ M, then we should have instead p| 

_ r (fp Vjl-K 1 

''^-7(2^^ (c.-i5;(p))V(r/2)2^^^7^- 

That is, we should allow the energy of the unstable particle to vary according to a Breit- 
Wigner function, rather than be fixed at = -E(p). (For simplicity of notation, we have 
left out a time-dilation factor of M/i?(p) that should multiply F; this will not affect our 
subsequent analysis, which is primarily concerned with the nonrelativistic regime.) If we 
take the limit F — 0, then the Breit-Wigner function becomes biuj — i?(p)), and we recover 
Eq. (|l]). On the other hand, if we take T ^ F ^ M, then the integral is dominated by the 
region near ~ T, and we have instead 

^ for T < F < M. (3) 



180 M2 

This is drastically different than the usual result, Eq. (|^); in particular, there is no expo- 
nential Boltzmann suppression. We see that this is essentially because the 93 particles that 
are being counted in Eq. (0) are far off shell, with energy near zero. 

While Eq. (^) may seem plausible, it leads to some surprising conclusions. Let us assume 
(following IHJ^) that the particle decays into two massless spin-zero x particles via an 
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interaction Tiint = ^fi(px^j thus we have T = /i^/327rM. Now suppose that we place a hot 
gas of light X particles in a large box, at a temperature T <^ M. The number density 
of X particles is = C(3)T^/7r^ and their energy density is = 7r^T^/30. After thermal 
equilibrium is established between ip and x particles, there should be a number density of 
ip particles given by Eq. (|^). The corresponding energy density is obtained by including 
an extra factor of E{p) in the integrand of Eq. (and not, as one might guess, an extra 
factor of uj). For T <^ M, this implies 

p^ = Mn^ ~ FTVM . (4) 

We now see that the ratio of cp energy density to x energy density is independent of tem- 
perature: 

p^/p^ ~ r/M . (5) 

This strikes us as odd, since there would not seem to be a source of virtual heavy ip particles 
in the limit of zero temperature for the massless x particles. 

The situation worsens for the case of stable, annihilating (p particles treated by MY [Q. 
For an interaction of the form Tii^t = ^A<y9^x^, MY find 

~ \{T/Mf'^T^ . (6) 

This implies ~ A(MT)^/^T^, and hence 

pjp^ ~ \{M/Tf'' . (7) 

Thus, for T ^ M, we see that the energy density in virtual heavy (p particles greatly 
exceeds the energy density in on-shell massless x particles. (This is not ruled out by energy 
conservation; the original temperature of the x gas would simply drop as the energy flows into 
virtual particles.) Eq. (|^ would seem to imply that (for example) the cosmic microwave 
background radiation is accompanied by a much larger energy density of virtual heavy 
particles. We believe that this is not a tenable proposition. 

Where, then, is the flaw in the MY analysis? Consider a system coupled to an environ- 
ment via an interaction, 

H = Hsys + Hcnv + -f^int , (8) 

where we assume that ifsys and Hem are positive semidefinite operators. We wish to deter- 
mine the energy of the system when it is in thermal equilibrium with the environment. The 
most obvious candidate for this energy is 

Esys = {Hsysjr , (9) 

where the angle brackets denote canonical thermal averaging with subtraction of the zero- 
point energy, 

Tr p-^/'^ 

(■■■)t= Tre-^/. -(01 •••10)- (10) 
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This definition of E^ys is the one used by JMY and MY. However, it is reasonable if and 
only if 

(i^int)T| < (i^sys)T . (H) 

If Eq. {\n\) does not hold, then the interaction between system and environment is effectively 
strong (no matter how small the coupling may be), and the appropriate division between 
system and environment is unclear. 

The analyses of JMY and MY are based on consideration of a Caldeira-Leggett model 
^ of coupled harmonic oscillators, grouped into terms according to Eq. (^. We will show 
below that in this model, at low temperature and weak coupling, 

(i^int)T ^ -2(i/sys)T • (12) 

We see that the negative interaction energy more than compensates for the system energy, 
which our qualitative arguments indicated was much too large. 

To demonstrate Eq. (0), we use the model presented in 0. A slightly different model 
was used in [Q; we have checked that Eq. (|12D holds in the model of as well. The model 

of i is 

H,y, = El c^c , (13) 

roc 

i/env = / dujujb^{uj)b{uj) , (14) 
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Hi^t= duo^Ja{u) c^b{uj)+b\u)c . (15) 

Here c and b{uj) are harmonic-oscillator operators with commutation relations [c, c^] = 1 and 
\p{u}),b\u}')] = 6{uj — u'), (j{uj) is a frequency- dependent coupling, and Uc is a lower cutoff; 
we assume ujc ^ Ei. The exact solution of this model involves changing to new variables 
B{uj) such that 

roo 

H = i/sys + i^cnv + ^^int = / duj iU B\uj)B{u:) , (16) 



where [B{uj), B'^{ijj'y\ = 5{uj — u'), and the original operators are given in terms of the new 
ones via 



c 



r du^l^)f{u)B{uj), (17) 

JljJc 

b{uj) = B{uj) + 0{a) . (18) 
Here the function f{uj) is given by 

uj — El + li(co') + t7Ta{ujj 

where 

U{uj) = P rdu'^j^ . (20) 
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The 0{a) term in the formula for b{uj) will not be needed; we will treat the coupling as 
weak, (j{uj) <^ Ei, and work to leading nontrivial order in a. This means we can neglect 
Il{uj) compared to Ei, and treat Ei as the renormalized single-particle energy; this point is 
thoroughly discussed in 

We now wish to compute (ifsys)T and (iJint)r- (We can also compute (ifenv)T, but the 
result is infinite, due to the infinite number of harmonic oscillators in the environment.) 
This is entirely straightforward; the formula we need is 



B^{uj')B{u;))^ = -^jj^ 6{uj' - u) . (21) 



Using Eqs. (Il||T3JT|j2T|) , we have 



{H,y,)T = El / duj' duj ^a{u')a{uj) r{u')f{u) {B\u')B{uj)), 
= E, r duja{uj)\f{uj)\' 



= L "^"^ (uj - E,Y + vrV2(.;) ^^V^ " ^^2) 

We see the similarity with Eq. (0). At high temperature and weak coupling, the region near 
uo El dominates, and we have 

(gsys)T ^ ^ fora(^)«T~Ei. (23) 

This is the same result that one would obtain for a noninteracting oscillator. On the other 
hand, at low temperature the low-tu region dominates, and we have 

(ifsys)T ^ du for T « a{u) « El. (24) 

We now turn our attention to the interaction energy. We begin by computing 



{c%{u))t= / du' ^a{uj') f\J){B\uj')B{uj))T + 0{a) 

= ^)nu) + 0{a) . (25) 

From here on we do not display the 0{a) correction. We then have 

roo I p -, 

{Hint)T= duj^a{uj) {c^b{uj))T + c.c. 

J UJc 

- r du 2(^-i^i)a(a;) 1 .261 

At high temperature and weak coupling, we get 

/•CO IrriLjj^ 1 

{H,^,)t^P duJ^-^ ,^ ^ ioia{uj)<€.T^Ei. (27) 

Jujc 1-1^ — 6 ' — i 



5 



This is smaller than (ifsys)T due to a suppression factor of (j{uj)/Ei. Thus the interaction 
energy is small compared to the system energy, as it should be. If, however, we consider low 
temperature and weak coupling, then we get 



Comparing with Eq. (^4]) gives us Eq. (|^). 

Clearly, then, the proper identification of the system energy becomes a key issue. We 
do not have a definitive resolution of this puzzle. However, we offer the following comments 
regarding the specific problem of WIMP relic densities. 

It seems to us that the WIMP relics that survive to the present day (and constitute the 
nonbaryonic dark matter) must be on-shell particles. If they are to be virtual, they must be 
produced via processes involving the present microwave-background photons; as discussed 
above, we find this to be untenable. Therefore, the correct computation is, we believe, 
of the number density of on-shell particles. This leads us back to the usual Boltzmann 
equation, which assumes all particles are on-shell. In fact, among the clearest derivations of 
this equation is the one presented by MY, before they go on to consider virtual effects. The 
usual Boltzmann equation of course leads to the standard result for the relic WIMP density, 
which we believe is correct. 
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